Quantum oscillations of spin current through a III-V semiconductor loop 
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We have investigated the transport of spin polarization through a classically chaotic semiconductor 
loop with a strong Rashba spin-orbit interaction. We found that if the escape time of a particle 
is long enough, the configuration averaged spin conductance oscillates strongly with the geometric 
spin phase. We predict a sizable rotation of spin polarization along its flowing path across the loop 
from the injector to the collector. We have also discovered a quantized universal spin relaxation in 
a 2D reservoir connected to such a semiconductor loop. 
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In the emerging field of spin electronics, the recent 
achievement of spin injection into paramagnetic semi- 
conductors jjj makes it an urgent task to control the 
spin current in semiconductor nanostructures. The spin 
current in a 2D channel of narrow gap III-V semicon- 
ductor can be manipulated by taking advantage of the 
strong spin-orbit splitting of the conduction electron en- 
ergy, because the mechanism of such splitting produces 
a spin precession which depends on the electron quasi- 
momentum. One example is the spin valve transistor , 
in which spin polarisation precesses in a 2D semiconduc- 
tor channel between a ferromagnet spin injector and a 
ferromagnetic spin collector. The measured resistance is 
determined by the angle of spin rotation along the prop- 
agating path. This angle can be varied by adjusting the 
spin-orbit interaction (SOI) strength in the semiconduc- 
tor with an external gate ^j. 

In this Letter we will investigate an interesting phe- 
nomenon which can be observed in a 2D semiconductor 
loop as shown in Fig. 1. It is well known Q that due to 
the SOI, when an electron travels along a closed path, its 
wave function accumulates an additional phase ip. If the 
SOI is a linear function of the electron quasimomentum, 
this phase depends only on the shape and the length of 
the path. In multiconnected conductors the effect of this 
phase on electron transport is similar to the Aharonov- 
Bohm (AB) effect. For example, in a disordered 2D ring, 
ip adds itself to the AB phase in the Aronov-Altshullcr- 
Spivak oscillation of the DC electric conductance B, as 
well as to the AB oscillation of the electric conductance 
mesoscopic fluctuations || . 

In this Letter, instead of the oscillation of the elec- 
tric current, we will study the effect of the spin phase 
on the quantum oscillation of the spin current. How- 
ever, we will assume the motion of electrons as ballistic 
along their classical trajectories, rather than a diffusive 
transport inside the loop [|. O ne such trajectory 
is schematically illustrated by the zigzag-lined path in 
Fig. 1, although it can be curved by a smooth random 
potential produced by the modulation doped impurities. 
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FIG. 1: A schematic plot of the loop sample. The zigzag line 
represents classical trajectory a, and the arrowed curve is the 
smooth Z-path. n(l) is a vector normal to the /-path. 



We assume that the motion of a particle inside the loop 
is classically chaotic, and so the quasiclassic approach of 
Ref. can be applied. Nevertheless, to calculate the 
spin current through the loop we need to generalize this 
method by taking into account the spin degree of freedom 
and the SOI. With this approach, we will calculate the 
average spin current. This implies that its mesoscopic 
fluctuations will be averaged out. The corresponding ex- 
perimental performance is, for example, to average the 
results measured with several gate voltage sweeps. We 
will also ignore the weak localization correction, which is 
small because our system has a large number of trans- 
port channels through the loop. We will show that the 
so calculated spin conductance oscillates as a function of 
the SOI strength, and consequently can be controlled by 
varying the gate voltage. We would like to emphasize 
that these quantum oscillations appear in the classical 
spin conductance, which is a drastically different phe- 
nomenon from the AB effect. The AB effect is abcent 
in the average spin current when the weak localization 
effects are ignored. 

Following the Landauer approach, to study the spin 
dependent conductance we define 
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where t"^ is the transmission amplitude of an electron 
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at Fermi energy Ep propagating from the channel n and 
the spin state a in the injector to the channel m and 
the spin state (3 in the collector, itself is the a(3- 
element of the matrix t mn which operates on spin states. 
The usual spin independent electric conductance is sim- 
ply g={e 2 /h) E n ,m,a,/3 l*nm| 2 - If a spin orinented along 
the x-axis is injected from the injector, and its orienta- 
tion becomes along the y-axis when the spin is collected 
at the collector, let g xy represent this spin current pass- 
ing through the loop. The matrix elements can then be 
written as 



where 
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where at are Pauli matrices with i=x, y, z. The spin ori- 
entation can be detected by measuring the polarization 
of the emitted light jlj]. In such an experiment, the po- 
larization matrix of the emitted photons can be derived 
if we know g^. 

In narrow gap III-V semiconductor quantum wells, the 
SOI is dominated by the Rashba interaction Q with the 
interaction Hamiltonian 



H x . 



a <x x p , 
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where p is the momentum operator, and the vector <x 
has components a x , a y , and a z . Then the quasiclassical 
expression of g a p 10 can be easily obtained as 

9aPl s^^J2 t o( a ) t o(b)SfS^\ (4) 
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where to(s) is the spin independent transmission ampli- 
tude for a classical trajectory labelled by s. One such tra- 
jectory is schematically plotteded in Fig. 1 as the zigzag 
line. The explicit expression of to(s) as well as the bound- 
ary conditions are given in Ref. 0. The spin evolution 
operator S a along the a-trajectory is defined as 



S a =T 
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where the integration is along the a-trajectory, m* is the 
effective mass, and z is a unit vector parallel to the z- 
axis. The symbol T means that the Pauli matrices must 
be ordered along the path. 

The g a f3-ys given in (||) must be averaged over the 
mesoscopic fluctuations, using the procedure described in 
Ref. ^. Such averaging may be considered as a tempera- 
ture effect, or, in accordance with the ergodic hypothesis, 
as an average over an ensemble of loops of slightly dif- 
ferent shapes. Because of the rapidly oscillating phase 
factors in the quasiclassical amplitudes to, after the aver- 
aging procedure, in (Q) only the terms with a=b remain. 
Accordingly, we obtain from (|^) the so averaged spin con- 
ductance {gij} as 
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Similarly, the so averaged electrical conductance is sim- 
ply (2e//i)^ a |io(a)| 2 j and is spin independent. 

The evolution matrix can be parametrized using its 
property that it is a SU(2) representation of 3D rota- 
tions. In fact, we can express S a as a time ordered prod- 
uct of infinitesimal rotations corresponding to small shifts 
dr along the trajectory a. Each infinitesimal rotation 
is along the axis drxz through an angle 2\dr xz\/L so . 
These infinitesimal rotations are represented by opera- 
tors exp[(— i/L so )(dr x z)<x], and they sum up to make 
S a for a finite rotation through the angle 2ip a around a 
unit vector N Q . Hence, the evolution matrix in can 
be represented as 



S a 
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We should notice that tp a and N a are uniquely deter- 
mined by the geometric shape and the length of the tra- 
jectory a. 

From now on we will consider a particular sample ge- 
ometry that the area occupied by the 2D electron gas in 
the loop is much less than L^ , and the linear dimension 
of the loop can be larger than L so , where L so =h/am* . 
In other word, both the upper path and the lower path 
of the loop are narrow. In this case one can show that 
each trajectory a in (^), as indicated by the zigzag line 
in Fig. I, can be replaced by a smooth trajectory, which 
is shown in Fig. I as the arrowed curve. We will label 
this smooth curve as /-path. Let 9 a be the area enclosed 
by the classical trajectory a making one turn around the 
loop. Then, the area enclosed by the Z-path is the average 
of 9 a over classical trajectories. Deviations of real paths 
from the Z-path can be treated perturbatively, which will 
be reported elsewhere. 

After the trajectory a is replaced by the Z-path, the 
evolution matrix (Q) becomes a simple function of the 
number w of windings the trajectory a makes around the 
loop until a particle escapes into the collector, w is posi- 
tive if the winding is counterclockwise. The correspond- 
ing evolution operator for the smooth Z-path, denoted as 
S(w), can be expressed as 



S(w) 
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Here 4>o and No are the 3D rotation parameters for the l- 
path in the lower half of the loop, and ip and N are the 3D 
rotation parameters for the Z-path around the complete 
loop. Using Eq. (||) we obtain the general dependence of 
the trace in (0) on the winding number 
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where the matrix elements M, 



(i) 



on the geometric shape of the Z-path. 
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Based on the above expressions, we can follow the ap- 
proach used in Ref. ^ to calculate the Aharonov-Bohm 
effect on mesoscopic electric conductance fluctuations in 
doubly connected classically chaotic loop. Let T be the 
time interval that a particle spends inside the loop, and 
To be such a duration for the shortest trajectory. Then, 
according to Ref. |^, we average the winding number with 
the Gaussian distribution function 
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where (3 is the system dependent dimensionless constant. 
For classically chaotic systems, it has been shown (71 that 
T obeys the distribution function V(T)—t^ 1 exp[— (T — 
Tq)/t], where r is the mean escape time of the parti- 
cles. We should point out that this kind of statistical 
approach is valid only for large r, such that the particle 
can travel around the loop many times before it escapes 
from the loop. Hence, we assume /c=^/2To/r/3<Cl. After 
averaging (10) over w and T we arrive at 
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A(tp) = K 2 /(K 2 + 4sin 2 ^). 
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The oscillation pattern of the spin current, caused by the 
A(ij)) term, gets sharper as k becomes smaller, and even- 
tually transforms into a periodic array of narrow peaks 
at positions ip=Trn. 

Eq. (Jl^) provides a general relation between the spin 
conductance and the spin phase ip. Here we will con- 
sider a specific example that the loop is nearly circular 
of radius cZ, and is symmetrically attached to two leads. 
For this geometry, by taking proper derivative with re- 
spect to the length of the Z-path, we obtain from (g) the 
differential equation 



(13) 



where the vector n(Z) is shown in Fig. 1. For a cicular 
Z-path, this equation is equivalent to the time-dependent 
Schrodinger equation for a | spin interacting with a uni- 
formly rotating planar magnetic field. The solution of 
this equation is well known JToj . Let us define r y=d/L so , 
and go the spin independent conductance of the loop. 
Referring to the coordinate system specified in Fig. 1, 
we obtain ^>=7ryl + 4 7 2 , and the nonzero components 
of the spin conductance as 



(9xx) 
(9xz) 
(9yy) 

(9zz) 



-go(n 2 /ip 2 ) [4 7 2 + A{ip) cos -0] 

- (ff»> = go^i: 2 H?) [1 - A{i>) cos 

— go A(tp) cost/' 

goty/^ 2 ) [1 + 4 7 2 A(V>) cos ip] . (14) 



From these expressions it is obvious that the oscillations 
of g can be observed as long as d>L so . From Ref. we 



evaluate L so ~3000 A in InGaAs/InAlAs quantum wells. 
Hence, the loop size must be larger than about 1 /j,m 
but less that the electron dephasing length, which can be 
very long at low temperatures. 

So far we have considered the spin transport between 
the injector and the collector through the loop. The 
mechanism of transport process involves spin diffusion 
driven by the difference of spin polarizations between 
the lead attached to the injector and the lead attached 
to the collector. Such a transport is represented by the 
spin conductance (Q) which, is determined by the spin 
dependent transmittance of the loop. Besides this pro- 
cess the SOI in the loop gives rise to a spin dependent 
reflectance. A particle which enters the loop from the 
injector-connected lead with a given spin orientation can 
be reflected back into the same lead with an opposite spin 
direction. This provides an additional relaxation process 
of the spin polarization in this lead, and the oscillatory 
dependence of such a relaxation on the spin phase ip is 
expected. A very suitable system for studying this relax- 
ation mechanism is a reservoir connected to a loop via a 
point contact which has M transmitting channels. The 
reservoir needs not to be very big. At a certain time, a 
nonequlibrium spin polarization E=(ATf— -ZVjJ/2 is created 
in the reservoir, where N a is the number of particles with 
spin projection a onto the quantization axis. Let be 
the reflectance associated to the electron spin flip reflec- 
tion. Then, the time rate of change of £ is given by 
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where \i„ is the chemical potential of the er-spin state, 
and t s is the spin relaxation time within the reservoir. 

Let the reservoir be a 2D degenerate electron gas of 
volume V, with the density of states at the Fermi level 
N(E F )=m* /2irh 2 . Then, (^-^=2^/ N(E F )V, and wc 
find from ( [l5| ) an additional spin relaxation rate in the 
reservoir 



r=(h/m*V) R U 



(16) 



due to its connection to the loop. Followinging the same 
quasiclassical approach which leads to the results (H), 
we obtain 



R 



= Af( \ (e™^* 



(17) 



where < • ■ ■ >„ is an average over the winding num- 
ber to. If the direction of spin polarization is parallel 
to N, then, the evolution matrix in ( O ) is diagonal, and 
hence T=0. On the other hand, if the spin polarization is 
perpendicular to N, we get |(exp[iw(/;N<T])|J 2 =sin 2 wtp. 
Consequently, averaging over w gives the final form for 
the relaxation rate 



F = (Mh/m*V) [1 - A(i/>)] . 



(18) 
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For small k the function oscillates between 

zero for -0=ri7r and a value very close to 1 for ?/>=(n+i)7r. 
Hence, around ip={n+^)it the spin relaxation rate T is 
determined mainly by the ratio Af/V. Taking a typ- 
ical value V—4 /im 2 , and the electron effective mass 
77i*/m=0.03 as in InAs, we obtain r~A/H0 9 sec -1 . The 
"quantum" of relaxation time, 10 -9 sec, is comparable 
to the intrinsic relaxation time r s . We would like to em- 
phasize that at this fixed spin phase ^=(n+^)iT, the re- 
laxation rate is a universal value in the sense that it does 
not depend on a shape and other parameters of the loop. 
On the other hand, we should also point out that (|l^) is 
valid only in the quasiclassic limit when A/"S>1. There- 
fore, one can not extrapolate this result to A/"~l where 
this quantization effect may be observed experimentally. 

In order to be able to observe the spin relaxation due 
to the SOI in the loop, it is necessary to have a long in- 
trinsic spin relaxation time in the reservoir. Therefore, 
the SOI in the reservoir must be sufficiently small while 
the SOI in the loop must be sufficiently large. Such a 
condition can be created by proper modulation doping, 
resulting in a highly asymmetric confinement potential 
near the loop but a nearly symmetric confinement po- 
tential in the vicinity of reservoir. A strong Rashba SOI 
will then appear only in the loop. 
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